Scale-free Network on Euclidean Space Optimized by Rewiring of Links 
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A Barabasi- Albert scale-free network is constructed whose nodes are the Poisson distributed 
random points within a unit square and links are the straight line connections among the nodes. 
The cost function, which is the total wiring length associated with a such a network defined on a 
two dimensional plane is optimized. The optimization process consists of random selection of a pair 
of links and rewiring them to reduce the total length of the pair but with the constraint that the 
degree as well as the out-degree and in-degree of each node are precisely maintained. The resulting 
optimized network has a small diameter as well as high clustering and the link length distribution 
has a stretched exponential tail. 

PACS numbers: 02.50.-r 89.20.-a 89.75.-k 89.75.Hc 



The Internet is a very complex network connecting the 
large number of computers around the world Q ■ The 
nodes of this network may be interpreted as the routers 
and links as the cables connecting computers. The net- 
work can also be described in the inter-domain level 
where each domain is represented by a single node and 
each link is an inter-domain connection. Such a network 
is well described by a graph consisting of a set of vertices 
and another set of edges among the vertices. Without 
assigning any weight with the links between the nodes 
only the topological structure of the Internet is mean- 
ingful. Study of Internet's topological structure may be 
important for designing efficient routing protocols and 
modeling Internet traffic. 

Internet is one of the large class of real-world net- 
works that exhibit small- world and/or scale-free prop- 
erties, e.g., social networks Q, biological networks 0,|5|, 
electronic communication networks 00 etc. Quantities 
that characterize a network of N nodes are the diam- 
eter V(N) which measures the topological extension of 
the network, the clustering co-efficient C measures the 
local correlations among the links of the network and 
the nodal degree distribution V(k). In a small- world 
network (SWN) §5], the diameter V(N) of the network 
scale logarithmically with N where as for a scale- free net- 
work (SFN) the degree distribution has a power law tail: 
V(k) oc fc~ 7 . Barabasi and Albert (BA) showed that a 
growing network with preferential link attachment prob- 
ability is a SFN with 7 = 30. 

Waxman first studied probabilistic graph models of In- 
ternet where links have weights which are their physical 
lengths 8] . The link length distribution in such networks 
decays exponentially: D(£) ~ exp(— l/i ). Faloutsos et. 
al. observed that the out-degree distribution of Internet 
follows a power law tail 1]. Yook et. al. observed that 
the distribution of routers of North America is a fractal 
set and the link length distribution is inversely propor- 
tional to the link lengths |fj|. It is suggested that in the 
growing Internet when a new node is becoming member 
of the network two competing factors control the deci- 
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FIG. 1: Two possible rewirings of a pair of links to reduce 
the total length of the two links. 



sion to which node of the already grown Internet the 
new node will be connected. The factors are the degree 
fej of the existing node i and in general the a-th power 
of the length £ of the link connecting the new node and 
the node i. The preferential attachment probability for 
the i-th node is therefore: Hi oc ki£ a . 

Recently it has been argued that such a network is 
scale-free for all values of a > a c = 1 — d in d dimen- 
sion and the link length distribution generally follows a 
power law D(£) ~ £ 5 where 5(a) = a + d - 1 pH EJ- 
For a < a c the degree distribution decays stretched ex- 
ponentially but D(£) still maintains a power law where 5 
saturates at — d — 1. The limit of a — > — oo is interesting 
where each node connects only to its nearest earlier node. 
In a regular network in the form of a linear chain sim- 
ilar studies have been done |l2{. An interplay between 
the preferential attachment and the link length selection 
within an interaction range for the Euclidean networks is 
studied in [l3j . 

In this paper we associate a cost function associated 
with such networks. Each link of the network has the 
cost equal to its Euclidean length £ and therefore the 
cost function of the whole network is the total length 
of all the links of the network. The question we ask is, 
how can one construct a small-world scale-free network 
with minimal cost? To study this we start generating a 
N node BA SFN on a two-dimensional plane. Links are 
then interchanged to reduce the cost function keeping the 
topology i.e. the degree value of each node intact. The 
optimization of the wiring length of networks on lattices 
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FIG. 2: Snapshots of the optimized network generated by the 
random rewiring process for iV=512 and m=2. (a) The initial 
network with £ « 395 (b) £ fa 304 at £=1000 (c) £ » 170 at 
i=10000 and (d) £ » 68 at t=10000000 where t is the number 
of rewiring trials. 



has been studied in |l4j . 

We start with constructing a Barabasi-Albert SFN 
embedded in the Euclidean space as follows. Let 
(xi, X2, xjf) and (y%, y2, Dn) be the independent 
identically and uniformly distributed random variables 
on the interval [0,1]. To construct one random configu- 
ration of the network let a specific set of values of the 
N pair variables (2:2, 2/2)) ( x NiUn)} be the 

co-ordinates of the set of N points on the unit square 
representing the set of nodes of the network with serial 
numbers i = 1 to N assigned to them. We use first m + 1 
points with serial numbers i = 1 to m + 1 to construct a 
m + 1-clique by connecting each point with rest of the m 
points. Then following the serial numbers new points are 
added to the network one after another and each node 
is connected to randomly selected m distinct previous 
nodes. The probability to link the new node with serial 
number j to a previous node i is linearly proportional to 
its degree, fcj. The network thus constructed up to N 
nodes is exactly the B A network Q • At the same time it 
is a small-world network, i.e. the diameter T>(N) of the 
network measured by the maximal distance between an 
arbitrary pair of nodes grow as log(iV). In this paper we 
restrict ourselves to m = 2. 

Let a denote the symmetric adjacency matrix of size 



N x N for our network such that ou 



1 if there is a 



link between the pair of nodes i and j and otherwise. 



FIG. 3: (a) The link length distribution D(i) in the optimized 
network: time-ordered (solid line) and random (dashed line) 
scaled by the average link length (£(N)) for iV=1024. (b) 
The average length (£(N)) varies as (N/lagN)^ in the opti- 
mized network: time-ordered (circles) and random (squares) 
rewirings. 



Let iij denote the shortest Euclidean distance between 
the pair of nodes i and j taking into account the periodic 
boundary condition. Therefore when = 1, £{j is the 
length of the connecting wire of the link between i and j. 
The total cost function C(N) is therefore the sum over 
all link lengths of the network i.e., C(N) = H^jaijiij. 

For the convenience of discussion wc define a degree 
vector similar to the contact vector generally used in the 
polymer physics. Our degree vector c describes the topo- 
logical connectivity of the network and has N elements 
Cj = hi, the degree of the i-th node. In the initial BA 
scale-free network one can associate a notion of time as 
if nodes are introduced one at each time unit. Therefore 
the links of the node i introduced at time i are divided 
into two groups 'outgoing' and 'incoming'. Each node 
has only k out = m outgoing links connected to m other 
nodes which are older than this node and it can be con- 
nected to k ln = k — m other nodes which are younger 
than this node. Consequently the degree vector can be 
split into two other degree vectors c out and c ln such that 

c out = k out and c m = k in and c = c out + c i„ 

Next, we perform the optimization dynamics to min- 
imize the total cost function C(N). The optimization 
dynamics conserves the number of links in the network, 
in addition it not only maintains the same degree vector 
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FIG. 4: The average diameter of the network T>(N) after op- 
timization, as a function of number of nodes N in the network: 
(a) for the time-ordered exchange V(N) — A + B log N where 
A=1.22 and .8=1.11 where as (b) For the random exchange 
V(N) ~ TV" where v = 0.31 ± 0.04. 

c but also c out and c m separately and thus ensures that 
the degree distribution of the network remains exactly 
same as it is before the optimization process starts. We 
call it as 'time-ordered' rewiring. One trial of rewiring in 
the optimization scheme consists of selecting four nodes 
ni,ri2,n3 and 77.4. The first node rii is randomly selected 
from the set of N nodes. 772 is selected randomly from the 
ki neighbours of n\. Similarly 713(7^ n\ =^ 77,2) is selected 
randomly from N nodes and 77,4(7^ n\ ^ 712) is again one 
of the &3 neighbours of 77.3. The move must maintain the 
conservation of link numbers as well as degree distribu- 
tion. We replace the link pair 77,177,2 and 77,3774 by another 
pair of links if either of the following two conditions is 
satisfied: 

i. if a%3 — ei24 = and £\2 + £34 > £13 + £24 we link 
77,177,3 and 77,2774. 

ii. if <2i4 = <223 = and £12 + £34 > £14 + £23 we link 
711774 and 77,277,3. 

If both are satisfied we accept one of them with probabil- 
ity 1/2. If only one is satisfied we accept that (Fig. 1). If 
none of the two is satisfied we go for a fresh trial. We also 
study a second type of rewiring process where only the 
total degree vector c is maintained but not individually 
c out and c ln . Here in the final optimized network a par- 
ticular node may have all neighbours which are younger 
than this node. We call this process as the 'random' 
rewiring method. In Fig. 2 we show how an initially 
complicated network becomes less messy with increasing 
number of rewiring trials. 

Since we accept the move only if the total rewired 
cost is reduced the trial is similar to the zero tempera- 
ture Monte Carlo dynamics. The total cost C(N) mono- 
tonically decreases with the number of successful trials 



and the number of un-successful trials between succes- 
sive accepted moves increases. We typically try around 
10(miV) 2 such trials so that the plot of C(N) with loga- 
rithm of the number of trials nearly reaches a plateau. 

From each point one can measure N — 1 distances and 
if these distances are sorted in an increasing sequence, 
one has the first neighbour distance, second neighbour 
distance, ... (N — l)-the neighbour distance etc. It is 
known that the average 77-th neighbour distance (R^) 
varies as N~ x / 2 if 77/iV << 1 and it is of the order of 1 
when n/N ~ 1 in the limit of N — > 00 [lfij . There is no 
other variation like N~ x when x is neither nor 1/2 but 
in between. 

In the optimized network the links are not necessarily 
a fixed (n) neighbour distances but a complex mixture 
of many neighbour distances. More elaborately it is ex- 
pected that many of the links of the optimized network 
are first neighbour distances, less numbers are second 
neighbour distances, less numbers are third neighbour 
distances etc. In the optimized network we first calcu- 
late the probability density of the link length distribution 
D(£). This distribution on scaling by the average link 
length (£(N)) is nearly the same for the time-ordered as 
well as the random rewiring processes. Contrary to the 
expectation this distribution has a maximum and it fits 
very well to a functional form D(£){£(N)) ~ x a e~ axf> 
with x — £/(£(N)). The fit on a linear scale gives 
a = 1.4,1.1 and (3 = 0.8,0.9 approximately for the 
time-ordered and random rewiring processes respectively. 
The network has Nt — 2N — m — 1 links and there- 
fore C(N) = N e (£(N)). We plot in Fig. 3 (£(N)) with 
N / log N and observe excellent straight lines on a double 
logarithmic scale. Therefore (£(N)) ~ (N/ log N)^ where 
/i = 0.46 and 0.52 with an error of 0.05 approximately 
for the time-ordered and random rewiring processes re- 
spectively. 

The topological size of the network is measured by the 
diameter of the network. The distance d\j between an 
arbitrary pair of nodes i and j is the number of links 
on the shortest path connecting the two nodes. The di- 
ameter d m is the maximal distance on a network. The 
average diameter T>(N) represents the configuration av- 
eraged maximal distance (d m ). Variations of the average 
diameter of the optimized network with the network size 
is shown in Fig. 4. For the the time-ordered exchange 
V{N) = A + BXogN with A w 1.22 and B w 1.11 
where as for the random exchange V{N) ~ N u with 
v = 0.31 ±0.04 where the error 0.04 is estimated by the 
largest difference of the local slopes between successive 
points and the mean slope. Although the degree distri- 
bution remains scale-free in the optimized networks gen- 
erated by both the time-ordered and in random rewiring 
procedures, the first network retains some long range con- 
nections due to the constraint that both c out and c ln 
are strictly maintained where as in the second network 
by random rewirings essentially all the connections are 
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FIG. 5: The clustering co-efficient C(k) as a function of the 
degree k shows a {k/{ln(k)} 1 ^ 2 }~ b behaviour with b ~ 0.94 
and 1.1 for the time-ordered (circles) and random (squares) 
rewiring processes respectively. 



of the Internet. While some links in Internet are the 
cable-less (microwave) links, many connections are made 
by real physical Ethernet cables. Therefore the question 
of optimizing cost of the total wiring length of the net- 
work arises naturally which is the main point of study in 
this paper. An optimized geographical embedding algo- 
rithm for scale-free networks was recently studied inde- 
pendently 0| • Unlike in [J^ , our time-ordered optimiza- 
tion produces a statistically non-homogeneous network 
and preserves a significant number of long-distance con- 
nections, permitting the network diameter to still scale 
as log(iV) as N — > oo. We also obtain a stretched- 
exponentially decaying tail of the link length distribution 
in the optimized network which is unlike the power-law 
tail observed by Yook et.al. || and closer to the Waxman 
result 0. 
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local, i.e. typically a node has all neighbours within a 
spatial distance of the order of ~ TV -1 / 2 . 

The local correlation among the links is measured by 
the clustering co-efficient. The clustering co-efficient d 
of the i-th node is measured by the ratio of the number of 
links e, within the ki neighbours of the i-th node and the 
number of links ki{ki — l)/2 if the ki nodes have formed 
an fci-clique i.e., Cj = 2ei/{ki(ki — 1)}. The clustering 
co-efficient of the whole network C(N) is (C). Also the 
average clustering co-efficient for the set of nodes of de- 
gree k is defined as C(k). In general both these clustering 
co-efficients may decrease as power laws: C(N) ~ 7V~ a 
and C(k) ~ k~ b . In our case we start from the initial BA 
network where it is known that a w 3/4 and 6 = 00. 
We also calculate these quantities in the final optimized 
state. The total clustering coefficient is found to be inde- 
pendent of N and therefore a=0 where as unlike a simple 
power law for C(k) we get a power law with logarithmic 
correction. In Fig. 5 we plot C(k) with fc/{in(fc)} 1 ' 2 
and observe straight lines on a double logarithmic scale 
implying the variation as: 



C(k) ~ {k/{ln{k)Y 12 } 



1/21-6 



(1) 



where b w 0.94 and 1.1 for the time-ordered and random 
rewiring processes respectively. We cannot rule out the 
possibility that b is actually 1 for both the processes. 
Many networks and models show 6=1 0, ■ 

To summarize, we have studied a cost optimized net- 
work which has three main features of the real- world net- 
works e.g., it is a small- world network, it is a scale-free 
network and also it exhibits high clustering properties as 
well. We studied this network on the two-dimensional 
Euclidean space which should be relevant in the context 
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